Hexagonal and Square Flux Line Lattices in CeCoIn 5 
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Using small-angle neutron scattering, we have imaged the magnetic flux line lattice (FLL) in the 
d-w&ve heavy-fermion superconductor CeCoIns. At low fields we find a hexagonal FLL. Around 0.6 
T this undergoes what is very likely a first-order transition to square symmetry, with the nearest 
neighbors oriented along the gap node directions. This orientation of the square FLL is consistent 
with theoretical predictions based on the d-wave order parameter symmetry. 

PACS numbers: 74.60.Ec, 74.70.Tx 



Recently, a whole new family of heavy fermion super- 
conductors has been discovered. It includes pressure- 
induced superconductivity in CeRhlnsEI, ambient pres- 
sure superconductivity in CeIrIn 5 B, and superconductiv- 
ity with the highest known T c = 2.3 K, for any heavy 
fermion at ambient pressure in CeCoInsB. The CeMIns 
(M = Rh, Ir, Co) family exhibits several similarities to 
other correlated electron superconductors such as high- 
T c cuprates and crystalline organic metals: their crys- 
tal structure consists of alternating units of_CeIn3 and 
MIn2 stacked sequentially along the c axisoO, the su- 
perconducting state borders on,-a magnetically ordered 
phase giving rise to competitiacJGJ or coexistence of mag- 
netism and superconductivitycnU. Finally there .is ev- 
idence, from thermal conductivity measurementsa and 
NMRa indicating that CeCoIns is a (i-wave superconduc- 
tor, with line nodes along the [110] and [110] directions 
(d x 2_ y 2). Theoretically, cZ-wave pairin g is | expej cted to 
stabilize a square flux line lattice (FLL)MllilOij, which 
was indeed recently reported in the high-T c superconduc- 
tor Lai. 83 Sro.i7Cu04 + 5 (LSCO)li-i However, in the lat- 
ter case with an orienta-tioa-jptated 45° with respect to 
theoretical predictionaljO'E3E3. In addition, it is worth 
pointing out that studies of the FLL symmetry in LSCO 
as well as in YBa2Cu307_a (YBCO) are susceptible to 
potential complications in interpretation due to the or- 
thorhombic crystal structure, which, leads to formaion of 
twin planes which can pin the FLDla. On the other hand, 
the crystal structure of CeCoIns is tetragonal which ex- 
cludes twinning, and this material may therefore turn out 
to be a better example of a "typical" <i-wave supercon- 
ductor. 

Here we report FLL imaging in CeCoIns, obtained by 
small-angle neutron scattering (SANS). The FLL under- 
goes what appears to be a first-order, field driven tran- 
sition from a hexagonal to a square FLL is observed. A 



square FLL has not previously been observed in a heavy 
fermion superconductor, and furthermore this is the first 
example of a square FLL in a d-wave superconductor ori- 
ented with the nearest neighbor directions parallel to the 
nodal directions of the gap. 

The SANS experiment was carried out at the Dll 
small-angle neutron scattering diffractometer at the In- 
stitut Laue-Langevin, Grenoble, France. Single crystals 
of CeCoIns were grown from an excess indium fluxtl, and 
had a T c = 2.3 K and B c2 (0) = 5.0 T parallel to the 
c-axis. The sample was composed of four single crys- 
tals with thicknesses t = 0.16 — 0.2 mm mounted side by 
side, each of which was individually aligned. The rather 
thin samples were necessary, due to the strong absorb- 
tion of low-energy neutrons by In. The total mass of the 
sample was 36 mg. Incident neutrons with wavelength 
A n = 0.6 nm and a wavelength spread AA n /A n = 10% 
were used, and the FLL diffraction pattern was collected 
by a 64 x 64 (1 cm 2 ) position sensitive detector. For all 
measurements, the sample was field cooled to 50 mK in 
a dilution refrigerator insert, placed in a cryostat with 
a superconducting magnet. Magnetic fields in the range 
0.3 to 2 T were applied parallel to the crystalline c-axis, 
and background subtraction was performed using mea- 
surements following a zero-field cooling. 

In Fig. [I] we show FLL diffraction patterns for applied 
fields of 0.3, 0.6 and 2 T. The images were constructed by 
summing a number of measurements at different angular 
positions, in order to satisfy the Bragg condition for the 
different peaks. A clear evolution of the FLL symmetry 
and orientation is evident. At the lowest field, twelve 
peaks are observed evenly distributed on a circle in re- 
ciprocal space as shown in Fig. |l](a). This corresponds 
to two hexagonal domains oriented along the [110] or the 
[110] directions. For a hexagonal FLL oriented with re- 
spect to an underlying square crystal symmetry, the exis- 
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FIG. 1: (color) FLL diffraction patterns for CeCoIn 5 obtained at 50 mK and applied fields of 0.3 T (a), 0.6 T (b) and 2 T (c), 
after subtraction of background measurement. In addition the data are smoothed by a 3 x 3 boxcar average, and the center of 
the image is masked off. The arrow in panel (a) shows the orientation of the crystalline axes. 



tence of two degenerate domain orientations having equal 
population is expected. As the field is increased to 0.6 
T, a primarily square FLL is found [Fig. 0(b)], with the 
majority of the scattered intensity concentrated in four 
irregularly shaped peaks. Again the FLL is oriented with 
the nearest neighbor direction along [110] and [110]. The 
square FLL remains stable up to the highest measured 
field of 2 T [Fig. 0(c)], where the diffraction pattern now 
shows four distinct Bragg peaks. The symmetry and ori- 
entation of the FLL is shown schematically in real space 
in Fig. 0, and compared to the symmetry of the super- 
conducting gap. While preserving the nearest neighbor 
direction along [110], the transition from hexagonal to 
square symmetry can in principle be continous. However, 
at intermediate fields this would result in diffraction pat- 
terns with contributions from 4 sheared hexagonal lattice 
orientations, analogous to what was previously observed 
in YBCOE3. Such a distortion of the hexagonal FLL was 
not observed, and the transition from square to hexag- 
onal symmetry is therefore most likely discontinous, i.e. 
of first order. On the other hand, we cannot exclude 
that the transition to square symmetry is preceeded by a 
weak rhombic distortion of the hexagonal FLL, similar to 
what was observed in the borocarbideslia. However, such 
a distortion would not alter the order of the transition 
from being first order. Finally, a co-existence of domains 
having respectively square and hexagonal symmetry is 
expected and usually observed in a narrow field range 
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FIG. 2: Symmetry and orientation of the hexagonal (a) and 
square (b) FLL. The orientation of the nodes in the super- 
conducting gap are shown in (c). 



around a first-order transitiorS. We expect this to be 
the reason for the slightly disordered square diffraction 
pattern seen in Fig. 0(b), and hence take the correspond- 
ing applied field of 0.6 T to be at or close to the transition 
field. The most likely first-order FLL symmetry transi- 
tion, and the orientation of the square FLL along the 
nodes of the superconducting gap, are the main results 
of this report, and will be addressed in further detail be- 
low. 

SANS FLL imaging in CeCoIns is complicated by two 
factors, illustrated by the diffraction patterns in Fig. 
and the 2 T rocking curve shown in Fig. 0. The most 
limiting factor is the long superconducting penetration 
depth in this material, which results in a very small 
field modulation and hence low scattered intensity. This 
inevitably leads to imperfect background subtractions, 
which is seen as negative scattered intensity in the diffrac- 
tion patterns in Fig. (dark blue regions). Fig. shows 
the intensity of a FLL Bragg reflection, as the cryostat 
is gradually tilted and rotated in such a way that the 
scattering vector cuts the Ewald sphere at a right angle. 
The integrated reflectivity is given byll3 
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FIG. 3: Rocking curve for the CeCoIns (1,0) FLL reflection 
at 50 mK and 2 T. The curve is a fit to a Gaussian. 
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where 7 = 1.91 is the neutron gyromagnetic ratio, t is 
the sample thickness, 4>q = h/2e = 2067 Tnm 2 is the 
flux quantum, and q — 2tv^/{B /<po) — 0.1955 nm -1 is 
the calculated scattering vector for a square FLL and a 
field of 2 T. This expectation and the measured value 
of q = 0.19 nm -1 agree within 3%. Theriiux line form 
factor h(q) for a square lattice is given bylia 



h(q) 



00 



(2ttA) 2 



-ttB/B c2 



where the exponential factor represents the so-called core 
correction. Fitting the rocking curve to a Gaussian and 
using the area under the curve as the integrated reflec- 
tivity, together with the upper critical field -B C 2(0) = 5.0 
T,0 we obtain A = 247 ± 10 nm. This falls inside the 
range of pkk*^-iAo = 190 — 281 nm reported in the 
litteraturetLa-EJEil. Calculating the coherence length from 
the upper critical field, £ = ^/(<j)o/2TTB c2 ) = 8.1 nm, 
we estimate the GL parameter k = A/£ « 30, making 
CeCoIns a strongly type-II superconductor. 

The second complicating factor is the narrow rocking 
curve, necessitating a very precise alignment in order to 
obtain scattering. The fit to the data in Fig. || yields 
a width of 0.21° FWHM comparable to the experimen- 
tal resolution estimated to be 0.15° FWHM. In principle 
the rocking curve width can be used to determine the 
longitudinal correlation length or straightness of the flux 
lines, but with the width being close to the experimental 
resolution this should rather be taken as a lower bound. 
We find Aq L < 0.21°(7r/180°)(7 = 7x 10~ 4 nm" 1 , and 
hence £l = 2/A^l > 3 /im. This is a large value corre- 
sponding to ~ 100 flux line spacings, and indicates very 
weak pinning in this material. 

We now return to the discussion of the symmetry and 
orientation of the FLL. In an ideal i s ot ropic type-II su- 
perconductor this will be hexagonalE3El However, if one 
evaluates the free energy difference between the hexago- 
nai-and square symmetry, this is found to be only about 
2%E3. A relatively weak anisotropy is therefore capa- 
ble of changing this delicate balance, leading to a dis- 
torted hexagonal or a square FLL. A number of theo- 
retical studies have addressed the effect of d-wave pair- 
ing on the structure and orientation of the FLL. As the 
field is increased or temperature decreased, they consis- 
tently find that a square FLL is stabilized, oriented with 
the neares t ne ighbor direction along the direction of the 
gap nodesE30't30. Determining the orientation of the 
hexagonal FLL is more difficult, since the energy differ- 
ence between the two configurations aligned 45° apart is 



very smalSEl. Ichioka et aZ.EH conclude that both the 
square and the hexagonal FLL are oriented with the near- 
est neighbors along the node direction, with a first-order 
transition separating the two symmetries. The transition 
field is predicted to be 0.15 x B c2 , which in the case of 
CeCoIns corresponds to 0.75 T at T — K. This is in 
agreement with our results concerning the orientation of 
the FLL as well as the nature of transition. Furthermore, 
their prediction of the transition field is in fair agreement 
with our estimate of ~ 0.6 T. 

In principle there is another mechanism that could be 
responsible for the FLL symmetry transition: A four-fold 
Fermi surface anisotropy combined with nonlocal electro- 
dynamics due to the finite coherence length. Theoreti- 
cally, this was studied extensively by Kogan et aL, who 
used nonlocal corrections to the London model to calcu- 
late the FLL free energy and thereby determine the-stablp 
configuration as a function of the flux line densityEirEilE3. 
This is the driving force behind the transition between 
a low-field (distorted) hexagonal and a high-field, square 
FLL seen in-t he rare-earth nickelborocarbidest§t3 as well 
as in VaSiEH. In the case of CeCoIn5 such an analysis 
has not yet been carried out. Band structure calcula- 
tions have been performed on the isostructural compound 
CelrlnsEl, and (partially) confirmed on CeCoIns hv mea- 
surements of de Haas - van Alphen oscillationsEa. The 
calculations show a Fermi surface with at least one sheet 
having a four-fold anisotropyO. However, this warps be- 
tween two orientations 45° apart, and at present it is not 
clear what implication this has on the FLL. 

To summarize, we have studied the symmetry and ori- 
entation of the magnetic flux line lattice in the <i-wave 
superconductor CeCoIns using small-angle neutron scat- 
tering. At low fields a hexagonal FLL was found, which 
undergoes which is most likely a first-order transition to 
square symmetry around 0.6 T. Though the possibility 
of a Fermi surface anisotropy combined with nonlocal 
effects can not be ruled out as the determining factor, 
our measurements agree well with the predictions for a 
pairing-symmetry driven transition. In particular, the 
nature of the transition, the field at which it occurs, and 
above all the orientation of the square FLL with the near- 
est neighbors alinged parallel to the node directions, are 
all consistent with being driven by the d-wave symmetry 
of the order parameter. 
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